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Comparison  of  Steady  State  Evaporation  Models 
for  Toxic  Chemical  Spills:  Development  of  a 

New  Evaporation  Model 


1.  INTRODUCTION 


The  United  States  Air  Force  handles  and  stores  a  number  of  toxic  and  hazardous 
chemicals.  Associated  with  this  activity  Is  the  threat  of  accidental  release  of  these  dangerous 
chemicals.  That  threat  applies  not  only  to  the  Immediate  area  of  the  spill,  but  to  locations 
downwind  of  the  spill.  To  determine  that  downwind  threat,  one  must  estimate  the  source 
strength  (evaporation  rate)  of  the  spilled  chemical.  One  spill  scenario  for  which  one  must  be 
prepared  is  that  of  a  liquid  chemical  spilled  onto  the  ground  so  that  It  forms  a  pool.  To 
estimate  the  evaporation  rate  from  the  pool,  data  Including  meteorological  information, 
properties  of  the  spilled  chemical,  characteristics  of  the  spill  site,  and  the  size  of  the  spill, 
must  be  readily  available.  The  procedure  for  estimating  source  strength  should  be  simple 
enough  to  run  on  a  microcomputer  with  a  minimum  of  knowledge  of  the  program  by  the  user. 

Several  evaporation  models  were  compared  by  Kunkel.1  The  Ille  and  Springer*  model  was 
identified  as  being  the  most  realistic  of  the  available  evaporation  models  because  It  allows  for 


(Received  for  Publication  28  November  1989) 

1  Kunkel,  B.A.  (1983)  A  Comparison  qf  Evaporative  Source  Strength  Models  fir  Toxic 
Chemical  Spills,  AFOL-TR- 83-0307,  ADA  1 39431 . 

a  Ille,  O.  and  Springer,  C.  (1978)  The  Evaporation  and  Dispersion  qf  Hydrazine  Propellants 
Jtom  Ground  Spills,  CEEDO-TR-78-30,  ADA  059407. 
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changes  In  the  pool  temperature  due  to  evaporation  and  solar  insolation.  Kunkel  modified  the 
Ille  and  Springer  model  so  that  it  includes  a  parameter,  n,  which  describes  the  wind  velocity 
profile.  This  parameter  was  referred  to  by  Kunkel  as  the  stability  index,  It  is  a  function  of 
atmospheric  stability  and  surface  roughness.  It  has  a  significant  impact  on  the  calculated 
source  strength.  A  100  percent  change  in  the  stability  index  number  generally  results  in  at 
least  a  43  percent  change  in  the  calculated  evaporation  rate.1,2  None  of  the  other  models 
described  in  the  comparison  by  Kunkel  considered  the  wind  velocity  profile.  The  Ille  and 
Springer  model  with  tills  change  will  be  referred  to  as  the  modified  Ille  and  Springer  model. 

Since  the  evaporation  model  comparison  by  Kunkel,  several  other  evaporation  models 
have  become  available  which  follow  the  same  general  calculation  procedure  as  the  Ille  and 
Springer  model.  This  procedure  is  a  steady  state  balance  of  all  sources  of  energy  that  add  to  or 
subtract  from  the  energy  of  a  pool  of  liquid  which  has  spilled  onto  the  ground.  Bach  of  these 
models  takes  a  different  approach  to  the  calculation  of  each  energy  input.  This  report  will 
serve  two  purposes:  (1)  It  will  compare  each  of  three  energy  balance  evaporation  models  and 
identify  the  most  appropriate  calculation  for  each  energy  input,  and  (2)  it  will  present  a  new 
energy  balance  evaporation  model  that  uses  the  most  appropriate  calculations  based  on  the 
comparisons. 

The  models  described  in  this  report  are  (1)  the  modified  Ille  and  Springer  evaporation 
model, 1>a  (2)  the  ADAM  model  liquid  pool  evaporation  source  calculation,3  (3)  the  Kawamura 
and  MacKay  evaporation  model,4  and  (4)  the  "New"  evaporation  model.  The  latter  is  a  new 
model  that  will  be  recommended  based  on  examination  and  evaluation  of  the  other  models. 

The  four  models  were  programmed  in  the  Basic  language  for  the  Zenith-248 
microcomputer.  The  ADAM  model  is  a  large  and  complex  model  including  many  types  of 
source  strength  calculations.  A  liquid  pool  evaporation  model  based  on  a  steady  state  energy 
balance  is  included  among  those  source  calculations.  The  ADAM  code  for  evaporation  of  a 
pool  of  liquid  was  extracted  from  the  ADAM  Fortran  code  and  translated  into  Basic  to 
function  as  a  stand-alone  model.  The  Kawamura  and  MacKay  model  was  taken  from  its 
Journal  article  description4  and  translated  into  a  Basic  code.  The  modified  Ille  and  Springer 
model  hod  already  been  written  in  Basic  code  for  the  Z-248. 

In  Section  2,  the  calculation  of  each  energy  term  from  the  modified  Ille  and  Springer, 
Kawamura  and  MacKay,  and  ADAM  models  are  evaluated,  Recommendations  are  made  for  the 
most  appropriate  calculations  for  a  new  evaporation  model  based  on  evaluations.  In  Section 
3,  the  new  model  is  summarized.  In  Section  4,  sensitivity  studies  are  presented  in  support  of 
the  new  model.  Concluding  remarks  are  made  in  Section  5. 


3  Raj,  P.K.  and  Morris,  J,A.  (1987)  Source  Characterization  and  Heavy  Gas  Dispersion  Models 
for  Reactive  Chemicals ,  AFGL-TR-88-0003  (I),  ADA  200121. 

4  Kawamura,  P.I.  and  MacKay,  D.  (1987)  The  evaporation  of  volatile  liquids,  Journal  qf 
Hazardous  Materials,  15:343-364. 
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2.  GENERAL  DESCRIPTION  OP  STEADY  STATE  ENERGY  BALANCE 


The  evaporation  rate  of  a  chemical  from  a  pool  surface  is  ultimately  a  function  of  pool 
surface  temperature.  The  equilibrium  surface  temperature  depends  on  many  avenues  of  heat 
transfer  to  and  from  the  pool.  These  include  solar  radiation  (Qtoi),  long  wave  radiation 
emitted  by  the  pool  (Qpol)  and  the  atmosphere  (Qatm).  convective  heat  transfer  from  the 
atmosphere  (Qh«)>  heat  conducted  from  the  ground  (Ogrd)<  and  heat  loss  due  to  evaporation  (Q«y), 
The  steady  state  temperature  is  the  temperature  at  which  the  sum  of  all  sources  of  heat  (energy) 
transported  Into  the  pool  exactly  balances  the  heat  transfer  out  of  the  pool;  that  is,  the  sum  of 
energy  terms  described  above  is  zero.  The  steady  state  energy  balance  Is  expressed  as; 


0«»1  +  0.tm  +  Opol +  ®he  +  °ev  +  ®  grd  "  0  total’ 


(1) 


where  Qtotni  ■  0  at  steady  state.  Many  of  these  energy  terms  can  be  expressed  as  a  function  of 
the  pool  surface  temperature.  The  equation  is  solved  Iteratively  for  the  pool  temperature, 
which  Is  then  used  to  calculate  evaporation  rate.  The  evaporation  rate  Is  proportional  to  the 
mass  transfer  coefficient  at  the  liquid  pool  -  atmosphere  interface  and  the  vapor  pressure  of 
the  chemical,  both  of  which  are  functions  of  the  pool  surface  temperature. 


2.1  Heat  Due  to  Solar  Radiation  (Ql0i) 

The  net  solar  radiation  reaching  the  spilled  liquid  depends  on  the  amount  of  cloud  cover, 
the  time  of  day,  and  the  geographical  location  of  the  pool,  The  last  two  parameters  are  used  to 
calculate  the  solar  altitude  angle. 

2.1.1  SOLAR  ALTITUDE  ANGLE  -  COMPARISONS 

In  the  modified  Ille  and  Springer  model,  the  solar  altitude  angle  is  input  by  the  user.  This 
amounts  to  a  visual  estimate  of  the  angle  of  the  sun  relative  to  the  surface  of  the  pool.  During 
cloudy  periods  this  could  be  difficult  to  do, 

The  ADAM  model  and  the  Kawamura  and  MacKay  model  calculate  the  solar  altitude  angle 
by  similar  methods.  They  both  start  with  the  following  equation  for  the  solar  altitude  angle, 
SA: 


sin  SA  ■  sin  LA  sin  D  +  cos  LA  cos  D  cos  SHA, 


(2) 
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where  LA  is  the  latitude,  O  is  the  solar  declination,  and  SHA  is  the  solar  hour  angle. 

The  models  differ  in  their  calculation  of  the  solar  declination  and  solar  hour  angle,  The 
ADAM  model  follows  the  procedure  of  Woolf,8  which  computes  the  exact  time  of  meridian 
passage  (true  solar  noon),  needed  to  calculate  the  solar  hour  angle,  The  Kawamura  and 
MacKay  model  uses  the  calculations  given  in  Lunde,6  which  sets  noon  equal  to  12,  The  solar 
declination  calculated  by  the  Kawamura  and  MacKay  model  is  likewise  a  simplified  version  of 
the  ADAM  model  solar  declination. 

2. 1 .2  SOLAR  ALTITUDE  ANOLE  -  RECOMMENDATION 

The  Kawamura  and  MacKay  and  the  ADAM  models'  solar  altitude  angle  calculations  are 
compared  in  Table  1 ,  The  differences  between  the  calculated  angles  for  the  two  models  are 
trivial.  Differences  in  solar  angle  as  small  as  these  are  not  expected  to  have  a  noticeable  effect 
on  the  calculated  evaporation  rate.  The  ADAM  model  solar  angle  calculation  is  slightly  more 
accurate  than  that  of  Kawamura  and  MacKay,  yet  it  requires  no  more  input  Information. 
Therefore,  the  ADAM  model  solar  altitude  angle  calculation  will  be  used  in  the  new  model, 


Table  1 .  Solar  Altitude  Angle  from  the  Kawamura  and 
MacKay  (KM)  Model  and  the  ADAM  Model, 


Time 


12:00 


12:00 


12:00 


12:00 


12:00 


12:00 


12:00 


14:00 


16:00 


18:00 


18:30 


19:00 


24:00 


KM  Solar  Angle 


8  Woolf,  H.M,  (1980)  On  the  Computation  of  Solar  Elevation  Angles  and  the  Determination  qf 
Sunrise  and  Sunset  Times,  National  Meteorological  Center,  Environmental  Science  Services 
Administration.  Hillcrest  Heights,  MO. 

n  Lund*-,  P.J.  (1980)  Solar  Thermal  Engineering,  John  Wiley  and  Sons,  New  York. 
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2.1,3  HEAT  FROM  NET  SOLAR  RADIATION  -  COMPARISONS 


In  the  ADAM  model,  the  net  heat  due  to  solar  radiation  Is  calculated  as  follows: 


Q#ol-RS(l-alb).Jm'V1, 


(3) 


where  alb  is  the  liquid  pool  albedo,  taken  to  be  0.14  (which  Is  a  standard  value  for  water  at 
Intermediate  lolar  angles),7  and  RS  is  the  net  radiation  per  unit  area  reaching  the  pool 
surface.  It  is  calculated  from  the  following:8*®10 


RS  « 


090  sin  (SA)-30^  _  ^  cp,3,4j 


(4) 


where  SA  is  the  solar  altitude  angle  and  RATI  is  the  indirect  radiation.  CT  is  the  cloud 
transmissivity  (related  to  thickness)  and  CF  is  the  fraction  of  sky  covered  by  clouds.  The 
diffuse  sky  radiation,  RATI,  is  a  function  of  solar  altitude  angle  and  is  calculated  from  the 
following: 


RATI  -  0.694  +  0.00349  SA  if  19.4  £  SA  <  42,  or 

RATI  m  0.49  +  0.014  SA  If  SA  <  19.4,  or 

RATI  «  0.84  lfSA*42. 


(8a) 

(Bb) 

(Be) 


In  the  modified  Ille  and  Springer  model,  the  net  heat  due  to  solar  radiation  for  a  cloudless 
day  is  calculated  from: 


7  MacKay,  D.  and  Matsugu,  R.S.  (1973)  Evaporation  rates  of  liquid  hydrocarbon  spills  on 
land  and  water,  The  Can.  J.  Ch.  Eng.,  81:434-439, 

8  Kunkel,  B.A.  (1988)  User's  Guide  for  the  Air  Force  Toxic  Chemical  Dispersion  Model 
IAFTOX),  AFGL-TR-88-0009,  ADA  199096. 

9  Holtslag,  A.A.M.  and  Van  Ulden,  A.P.  (1983)  A  simple  scheme  for  daytime  estimates  of  the 
surface  fluxes  from  routine  weather  data,  J.  Climate  AppL  MeteoroL,  22:517-529. 

10  Hasten,  F.  and  Czeplak,  G.  (1980)  Solar  and  terrestrial  radiation  dependent  on  the 
amount  ;\nd  type  of  cloud,  Solar  Energy,  24:177-189. 
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Q. 


sol 


1.1^4.x  .1?.-  (o.9c,cSA)  (sin  sa),  cal  m"ahr_1. 
RATI  v  '  v  7 


(6) 


In  this  equation,  1,164  x  10e  (cal  m-2  hr-1)  represents  the  solar  constant  and  0.9  represents  the 
transmission  coefficient  of  the  atmosphere  under  clear  conditions,  The  diffuse  sky  radiation, 
RATI,  is  the  same  as  In  the  ADAM  model.  If  clouds  are  present,  the  type  of  cloud  cover  is  input 
by  the  user  (one  of  8  choices  ranging  from  cirrus  to  fog).  Each  cloud  type  is  associated  with  a 
ratio  of  insolation  with  an  overcast  sky  to  Insolation  with  a  cloudless  sky.  Eq.  (6)  is 
multiplied  by  this  ratio  to  get  the  radiation  reaching  the  ground  when  clouds  are  present. 

There  is  no  provision  for  taking  amount  of  cloud  cover  into  account. 

In  the  Kawamura  and  MacKay  model,  the  net  heat  calculated  due  to  solar  radiation  is: 

Q#0l  -  4000  ( 1  -  0.0071  Ca)  (  sin  SA  -  0. 1 ),  kJ  m“3hr"  \  (7) 


where  C  is  the  cloud  cover  fraction  in  tenths,  and  SA  is  the  solar  altitude  angle.  Note  that  this 
model  does  not  account  for  the  cloud  thickness  (transmissivity), 

2. 1 .4  HEAT  FROM  NET  SOLAR  RADIATION  -  RECOMMENDATION 

The  solar  heating  equation  of  Kawamura  and  MacKay  (Eq.  (7)]  is  not  valid  for  solar  angles 
(SA)  less  than  5.74  degrees  (sin  SA  <  0.1).  Qao|  is  set  equal  to  zero  for  solar  angles  in  this 
range,  The  solar  heating  equation  of  the  modified  Ille  and  Springer  model  [Eq.  (6))  is  not  valid 
for  SA  ■  0.  In  those  cases,  SA  is  set  equal  to  1.  One  would  expect  a  aero  heat  input  due  to  solar 
radiation  during  the  night,  as  calculated  by  the  ADAM  model  and  the  Kawamura  and  MacKay 
model.  The  Kawamura  and  MacKay  model  does  not  consider  cloud  thickness;  the  modified  Ille 
and  Springer  model  does  not  consider  the  amount  of  cloud  cover.  The  cloud  type 
(thlckness/transmlsslvity),  as  well  as  cloud  cover,  are  expected  to  influence  the  solar  heating. 

Based  on  the  above  observations,  the  new  model  will  use  the  ADAM  model  method  of 
calculating  the  heat  coming  Into  the  pool  due  to  solar  radiation.  This  method  takes  into 
account  both  cloud  cover  amount  and  cloud  type  (thickness),  It  calculates  a  value  of  zero  for 
heat  due  to  solar  radiation  at  nighttime,  and  a  positive  value  for  solar  altitude  angles  greater 
than  zero. 
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a. a  Long  Wave  Radiation  Emitted  by  the  Atmosphere  (Q*tm) 


2.2.1  COMPARISONS 


The  total  amount  of  energy  radiated  by  a  black  body  per  unit  area  and  time  ia  given  by  the 
Stefan-Boltzmann  fourth  power  law: 

E-oT4  (8) 

where  T  ie  the  absolute  temperature  and  o  la  the  Stefan-Boltzmann  constant, 

8.87  x  10“*  J  m"a  a-1  K"1.  All  three  of  the  models  calculate  the  radiation  emitted  by  the 
atmosphere  using  the  following  equation: 


(9) 


The  additional  factor  et  represents  the  emisslvity  of  the  atmosphere  (■  1  for  a  perfect  black 
body),  e,  is  a  function  of  the  water  vapor  pressure  of  the  atmosphere. 

The  modified  Ille  and  Springer  and  ADAM  models  set  ea  equal  to  0.75.  The  Kawamura  and 
MacKay  model  calculates  atmosphere  emisslvity  via: 


et*(l  ~r)B, 


(10) 


where  r  is  the  reflectivity  of  the  pool  surface  and  is  set  equal  to  0.03.  B,  the  atmospheric 
radiation  factor,  is  a  function  of  cloud  cover  and  vapor  pressure  of  water  in  air,11  and  is 
determined  graphically.  A  graphical  determination  is  not  practical  for  a  computer  program. 
In  an  example  given  by  Kawamura  and  MacKay,  the  atmospheric  radiation  factor  is  equal  to 
0.84  at  18  mbar  water  vapor  pressure  and  20  percent  cloud  cover.  It  follows  for  this  example 
that  e*  •  0.81,  compared  with  0.75  assumed  for  the  modified  Ille  and  Springer  and  the  ADAM 
models. 


2.2.2  RECOMMENDATION 

The  calculated  evaporation  rate  is  not  expected  to  change  significantly  over  the  range  of 
typical  atmosphere  emlsslvities.  The  inconvenience  of  calculating  the  water  vapor  pressure 
and  determining  the  atmospheric  radiance  factor  from  that,  as  required  for  the  Kawamura  and 


11  Raphael,  J.M.  (1962)  Prediction  of  temperature  in  rivers  and  reservoirs,  Proc,  Amer, 
Soc,  C(c.  Eng.  J,  Power  Dtv.,  88:P02. 
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MacKay  model,  le  not  worth  the  minimal  benefit.  Therefore,  It  !•  reasonable  to  assign  a 
constant  value  of  0.75  to  the  atmosphere  emlsslvlty,  as  in  the  ADAM  and  the  modified  Hie  and 
Springer  models. 

2.3  Long  Wave  Radiation  Emitted  by  the  Pool  (Spot) 

2.3.1  COMPARISONS 

All  three  of  the  models  calculate  the  radiation  from  the  pool  from  the  Stefan-Boltsmann 
Law: 

Opd-V’^  <»> 


The  emlsslvlty  of  the  pool  (ep)  is  set  equal  to  0.95  in  the  modified  file  and  Springer  and  the 
ADAM  models,  and  0.97  in  the  Kawamura  and  MacKay  model.  These  values  are  approximately 
the  emlsslvlty  of  water.12 

2.3.2  RECOMMENDATION 

According  to  McAdams, 13  the  emlsslvlty  of  water  for  long  wave  radiation  ranges  from  0.95 
to  0,963  (for  temperatures  ranging  from  0  to  100  degrees  Centigrade).  For  the  new  model,  a 
pool  of  liquid  with  emissive  properties  similar  to  water  near  the  freezing  point  will  be 
assumed.  Thin  corresponds  to  an  emlsslvlty  of  0.95. 


2.4  Heat  Conducted  from  the  Ground  (figrd) 

2.4.1  COMPARISONS 

Heat  Is  transferred  from  the  ground  to  the  liquid  at  the  ground  surface  when  the  ground  is 
wanner  than  the  liquid.  This  heat  transfer  will  decrease  with  time  as  the  ground  temperature 
approaches  that  of  the  liquid  pool  temperature.  There  will  also  be  thermal  resistance  within 
the  pool  of  liquid  in  transferring  heat  from  the  bottom  of  the  liquid  layer  to  the  top,  where 
evaporation  occurs.  Heat  transfer  from  the  ground  to  the  liquid  is  driven  by  the  difference  in 
temperature  between  the  ground  and  the  pool: 


13  McAdams,  W.H.  (1954)  Heat  Tranamiaalon,  3rd  edn,,  McGraw-Hill  Book  Company,  Inc., 
New  York. 
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°«n|-h«<VTP> 


(12) 


The  modified  file  and  Springer  model  disregards  the  heat  transfer  from  the  ground.  It 
describes  only  the  heat  transfer  within  the  liquid  layer.  They  calculate  the  heat  transfer 
coefficient  from  the  equation: 


(13) 


whew  ki  *  thermal  conductivity  of  the  liquid  (cal  m^hr^H"1), 


Or  ■  Qrashof  number  (a  function  of  liquid  properties), 

Pr  ■  Prandtl  number  (a  function  of  liquid  properties), 

L  -  depth  of  the  liquid  pool  (m),  and 

A,  B,  C  -  constants  that  depend  on  the  heat  transfer  mode  (amount  of  turbulence), 


The  ADAM  model  assumes  that  heat  transfer  resistance  within  the  liquid  layer  is 
negligible.  It  describes  the  heat  transfer  coefficient  through  the  ground  as  follows: 


("  “« l) 


l/a 


.  J  m“as_1K' 


(14) 


where  kg  ■ 

«g  - 

t  - 

Cp,g  ■ 

P«  ■ 


thermal  conductivity  of  the  ground  y  iti"1  r1  K-1), 
thermal  dlffuslvlty  of  the  ground  fkg  Cp,g  pg)  (m3  s"1), 
time  after  spill  (s), 

heat  capacity  of  the  ground  (J  kg*1  K_l),  and 
density  of  the  ground  (kg  nr3). 


The  ADAM  model  authors  assert  that  the  natural  convective  effects  within  the  thin  liquid 
layer  are  negligible  compared  with  the  heat  transfer  through  the  ground. 
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The  Kawamura  and  MacKay  model  indudea  both  conduction  of  heat  from  the  ground  and 
heat  transfer  through  the  liquid  layer  in  an  overall  ground-liquid  heat  transfer  coefficient, 
Ugrd<  It  is  calculated  using  an  electrical  resistance  analogy: 


(18) 


hg  Is  the  coefficient  of  heat  conduction  through  the  ground,  and  hi  is  the  liquid  heat  transfer 
coefficient. 

Kawamura  and  MacKay  calculate  a  time-dependent  heat  transfer  coefficient  as  follows: 


(18) 


where  kg,  otg,  and  t  are  as  defined  for  the  ADAM  model. 

The  Kawamura  and  MacKay  liquid  heat  transfer  coefficient  is  estimated  from: 


h 


1 


(17) 


where  d  is  the  average  depth  of  the  liquid  pool,  approximately  equal  to  the  initial  pool  depth 
divided  by  2.  ki  is  the  liquid  thermal  conductivity  of  the  liquid  chemical  (kJ  nr1  -  hr*1  -  K'1), 
expressed  as  a  function  of  pool  temperature.  $  is  a  liquid  resistance  factor  for  which  they 
derive  the  following  empirical  equation: 


♦  ■  1/{1  +  exp [-0.06  (TB- 70)]}, 


(18) 


where  TB  is  the  liquid  chemical  boiling  point  in  degrees  C. 
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2.4.2  RECOMMENDATION 

In  the  Kawamura  and  MacKay  model,  hg  is  a  factor  of  2  higher  than  the  ADAM  model  hg. 
The  reference13  cited  by  Kawamura  and  MacKay,  uses  the  following  equation  for  the  transient 
heat  flow  in  a  semi-inflnlte  solid: 


% 


.-kA(V,T.) 

("Vr 


(19) 


where  qo  ■  heat  transfer  rate  at  the  solid  surface,  J  s"1, 

k  ■  thermal  conductivity  of  the  ground,  kj  nr1  -  s"1  -  C_l, 

A  ■  area,  ma, 

ag  -  k  p  C  ■  thermal  diffusivity  of  the  ground,  ma  s“l, 
p  ■  density,  kg  m*3, 

C  ■  speclflc  heat  capacity,  kJ  kg'1  -  C~l, 
t  ■  time  constant,  s, 

T0  ■  surface  temperature  of  the  ground,  and 
Ti  ■  initial  temperature  of  the  bulk  ground. 

A  number  of  other  references  provide  a  similar  analysis  of  the  time-dependent  heat  flow  in  a 
semi-infinite  solid.14'15  Therefore,  the  new  evaporation  model  will  describe  the  time- 
dependent  heat  transfer  coefficient  as  derived  In  the  aforementioned  references.  This  is  the 
same  calculation  as  In  the  ADAM  model. 

Both  the  ADAM  and  Kawamura  and  MacKay  models  assume  constant  values  for  the  ground 
properties.  But  those  properties  can  vary  substantially  depending  on  the  type  of  soil  or  surface 
on  which  the  chemical  has  spilled.  Table  2  lists  various  ground  types  and  their  corresponding 
thermal  conductivities  and  thermal  diffuslvltles.  In  the  new  model,  an  option  will  be  added 
whereby  the  user  can  choose  the  type  of  soil  or  surface  on  which  the  chemical  has  spilled. 


13  Holman,  J,P.  (1976)  Heat  Transfer,  4th  edn.,  McGraw-Hill  Book  Company,  Inc.,  New  York. 

14  Bird,  R.B.,  Stewart,  W.E.,  and  Lightfoot,  E.N.  (1960)  Transport  Phenomena,  John  Wiley  and 
Sons,  Inc.,  New  York. 

15  Schneider,  P,J.  (1965)  Conduction  Heat  Tranajer,  Addlson-Wesley  Publishing  Co., 
Cambridge,  MA. 
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Table  2.  Thermal  Properties  of  Different  Ground  Types. 


I _ I 


Soil/Surface 

Type 

Thermal 
Conductivity 
Wnr1  K"1 

Thermal 

Conductivity 

Reference 

Diffusivlty 

m*/s 

Diffusivlty 

Reference 

0.138 

16 

1.1  x  10“T 

• 

Wet  Soil 

0.687 

17 

1.9  xiO"7 

17 

Dry  Sand 

0.326 

17 

2.74  xUT® 

17 

Wet  Sand 

1.128 

17 

4.92  x  10“« 

17 

Concrete 

1.279 

17 

4.92  x  10“« 

17 

‘Estimated  from  a  ■  k/(p  Cp)  for  soil  organic  matter. 


As  shown  in  the  example  below,  the  thermal  resistance  of  the  liquid  as  calculated  by  the 
Kawamura  and  MacKay  method  is  non* negligible  compared  with  the  ground  heat  conduction, 


Example 

Let  chemical  «  N3O4  (TB  -21.2,  C) 
t  ■  10  min  -  600  s 
kg  -  0.68?  W  nr1  K“»  (wet  sod) 
oig  ■  1.9  x  10“r  (wet  soil) 

hg  -  0.687  /  (x  x  1 .0  x  lo-7  x  600)1'*  ■  34.7  W  nra  K-» 
d  -  0,8  cm  a  0.008  in  (average  pool  depth) 

«  -  1/1 1  +  exp  (-0.06  (21.2  -  70)))  •  0,08076 
k!  a  0.132  Wm-‘  K-» 

therefore,  hi  ■  820  W  nr3  K'1  and  Ugnj  a  32.6  W  nr*  K“l. 


The  overall  ground/liquid  heat  conduction  is  7  percent  lower  than  the  heat  conduction  from 
the  ground  only.  For  a  deeper  pool,  the  reduction  in  the  ground  heat  conduction  by  the  liquid 
is  even  greater.  For  Instance,  a  pool  whose  initial  depth  is  equal  to  8  cm  has  a  liquid  thermal 
resistance  that  lowers  the  overall  ground/liquid  heat  conduction  by  25  percent,  It  is  expected 
that  a  change  in  ground  heat  conduction  of  this  magnitude  can  have  some  effect  on  the 
calculated  evaporation  rate, 


16  Dean,  J.A,  (1988)  ed.  Lange's  Handbook  qf  Chemistry,  13th  edn.,  McGraw-Hill  Book 
Company,  New  York. 

17  Kothandaraman,  C.P,  and  Subramanvan,  S,  (1978)  Heat  and  Mass  Transfer  Data  Book, 
2nd  edn.,  John  Wiley  and  Sons,  New  York. 
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The  overall  procedure  euggeeted  by  Kawamura  and  MacKay  for  calculating  the  heat 
conducted  through  the  ground  and  liquid  pool  will  be  ueed  in  the  new  model,  with  the 
elimination  of  the  factor  of  2  in  Eq.  (16).  This  procedure  allows  pool  depth  to  be  selected 
(rather  than  assigned  a  constant  value,  as  in  the  ADAM  model)  and  takes  into  account  heat 
transferred  through  both  the  ground  and  the  liquid.  Eq.  (12)  then  becomes: 


°H-ur«CVTp) 


Jm“i 


-a.-i 


(30) 


2.0  Heat  Lots  Due  to  Evaporation  (Q«r) 

2.6. 1  COMPARISONS  -  GENERAL 

‘  *  ^ 

Vapor  transfe:  Lorn  the  surface  of  the  pool  to  the  surrounding  atmosphere  is  governed  by 
the  mass  transfer  coefficient  and  the  concentration  gradient  between  the  vapor  at  its  source 
and  in  the  ambient  surroundings: 

B*km(V0.)'k*m'V‘'  ,  01) 


E  is  the  mass  evaporation  rate  per  unit  area,  km  is  the  mass  transfer  coefficient  (m/s)  and  Cp 
and  C.  are  the  vapor  concentrations  at  the  pool  surfkce/air  interface  and  in  the  surrounding 
air,  respectively.  Assuming  C»  -  0  and  the  vapor  behaves  as  an  ideal  gas,  one  can  write 


kMP  . 
i  m  .  itt , 

RT_ 


kgm 


(32) 


where  M  (kg/mol)  is  the  vapor  molecular  weight,  PBtt  is  its  saturation  vapor  pressure  (Pa),  R  is 
the  ideal  gas  constant  (J  mol"1  K_1).  and  Tp  is  the  absolute  temperature  of  the  pool.  The  rate  of 
heat  loss  due  to  evaporation  per  unit  area  is  the  product  of  the  evaporation  rate  and  the  heat  of 
vaporisation,  Hv  (J/kg): 


(23) 
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2.5.2  MASS  TRANSFER  COEFFICIENT 

All  the  models  discussed  here  use  this  same  general  procedure  for  calculating  the 
evaporative  cooling  of  a  pool  of  liquid,  They  differ  in  the  calculation  of  the  mass  transfer 
coefficient.  Mass  transfer  is  a  fundamental  property  Important  to  the  evaporation  process.  It 
Is  therefore  Important  to  understand  the  details  of  the  calculation  of  this  property. 

In  the  modified  llle  and  Springer  model,  the  mass  transfer  coefficient  is  calculated  in 
terms  of  the  wind  speed,  U  (m/h),  and  pool  diameter,  X  (m)  (or  more  precisely,  the  distance 
across  the  pool  in  the  downwind  direction),  as  given  by  Sutton;18 


a-n  -n 

k  .  0.0292  Sc"0687  u*  +  n  X»  +  n,  m/hr. 

ttt 


(24) 


Sc  is  the  Schmidt  number,  a  dimensionless  number  calculated  by  dividing  the  kinematic 
viscosity  of  the  air/chemical  vapor  mixture  by  the  dlflxisivity  of  the  chemical  vapor  in  the  air. 
Mass  diffusivity,  like  mass  transfer,  is  calculated  differently  by  each  evaporation  model.  This 
will  be  discussed  later  in  this  section.  Hie  constant  n  is  chosen  to  describe  the  wind  velocity 
profile.  It  is  a  function  of  the  terrain  roughness  and  the  atmospheric  stability  (solar  altitude 
angle,  cloud  cover,  and  wind  speed)  and  is  a  very  sensitive  parameter,  n  has  been  estimated  by 
Sinedman-Hogstrom  and  Hogstrom.19  Ille  and  Springer  show  that  a  100  percent  change  in  n 
results  in  a  43  percent  or  greater  change  In  the  evaporation  rate. 

The  Kawamura  and  MacKay  model  calculates  mass  transfer  coefficient  in  basically  the 
same  way  as  the  modified  llle  and  Springer  model,  but  with  a  few  simplifications.  Rather 
than  selecting  n  based  on  terrain  type  and  atmospheric  stability,  the  model  chooses  a  constant 
value  of  0.20  for  n  to  represent  typical  atmospheric  and  terrain  conditions.  This  is  the  same 
calculation  as  the  original  llle  and  Springer  model.  Mass  transfer  coefficient  then  becomes: 


m 


0.0292  Sc“°e67Ua78X-011. 


(25) 


18  Sutton,  0,0,  (1903)  Micrometeorology,  McGraw-Hill  Book  Company,  Inc,,  New  York. 

18  Smedman-Hogstrom,  A.S,  and  Hogstrom,  U.  (1978)  A  practical  method  for  determining 
wind  frequency  distribution  for  the  lowest  200  m  from  routine  meteorological  data,  J,  qfAppL 
Meteorol,  17:942-954, 
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The  ADAM  model  calculates  the  mass  transfer  coefficient  as: 


k 


m 


ILbl  AQSL  (Re08 
X  cO^V 


15,50o),  m/s, 


(26) 


where  nm  is  the  kinematic  viscosity  of  the  pure  chemical  vapor,  X  is  the  pool  diameter,  Sc  is 
the  Schmidt  number  (as  defined  previously),  and  Re  is  the  Reynolds  number.  Re  is  calculated 
from  the  wind  speed,  pool  width,  and  viscosity  by: 


Re-U*. 


(27) 


The  ADAM  model  uses  the  kinematic  viscosity  of  the  pure  chemical  vapor  to  calculate  the 
Schmidt  number,  The  Kawamura  and  MacKay  model  uses  the  kinematic  viscosity  of  the  air, 
and  the  modified  Ille  and  Springer  model  uses  the  kinematic  viscosity  of  the  chemical 
vapor/alr  mixture.  The  Kawamura  and  Mackay  model  also  simplifies  the  Schmidt  number  by 
setting  it  equal  to  the  kinematic  viscosity  of  the  air  (rather  than  the  chemlcal/alr  mixture) 
divided  by  the  diffuslvlty  of  the  chemical  in  air. 

2.8.3  MASS  DIFFUSIVITY 

As  with  the  overall  mass  transfer,  the  models  differ  in  their  approach  to  the  mass 
dtffusivlty.  The  modified  llle  and  Springer  model  uses  an  estimation  method  for  dlfhisivity 
based  on  theory: 


DV 


•b 


0.00266  T3^  a , 

- ,  cm  /s, 


(28) 


where 


DVftb  ■  molecular  diffusion  coefficient  of  chemical  a  in  chemical  b, 


M 


tb 


“b 


(29) 
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M^.Mb  -  molecular  weights  of  a  and  b, 


o  -  characteristic  length  (A  )"(°,+ 
nD  -  diffusion  collision  integral. 

The  characteristic  length  and  collision  integrals  are  calculated  from  the  following; 

c  ■  1,18  V*/3  (Vfa  ■  liquid  volume  at  boiling  point),  (30) 

ftD  -  10  (“a43 lo«  Crp)/a<303  +  OMBOK).  (31) 


OMEOK  is  the  collision  integral  constant,  unique  to  each  chemical  and  not  widely  available. 

Tp  is  the  temperature  of  the  pool  of  liquid. 

Kawamura  and  Mackay  used  data  from  Holman18  for  the  moleoular  diffusion  coefficient  in 
the  example  that  they  gave.  This  value  was  given  for  a  specific  temperature  as  opposed  to 
calculating  it  for  each  temperature  iteration. 

The  ADAM  model  uses  the  method  of  Fuller  et  ala031,aa  to  calculate  the  molecular  diffusion 
coefficient: 


DV. 


0.00143  T 


1.78 


»b 


(33) 


30  Fuller,  E.N.  and  Qiddings,  J.C,  (1063)  A  comparison  of  methods  for  predicting  gaseous 
diffusion  coefficients,  J.  Gas  Chromatogr.,  3:222-227. 

31  Fuller,  E.N.,  Ensley,  K.,  and  Qiddings,  J  C.  (1969)  Diffusion  of  halogenated  hydrocarbons 
in  helium;  The  effect  of  structure  on  collision  cross  sections,  J,  Phys,  C hem.,  73:3679-3686. 

aa  Fuller.  E.N.,  Schlettler,  P.D.,  and  Qiddings,  J.C.  (1960)  A  new  method  for  prediction  of 
binary  gas-phase  diffusion  coefficients,  Ind,  Eng.  Chtm.,  68:18-27. 


16 


where  the  Ey's  are  the  molecular  diffusion  volumes.  The  other  parameters  have  been 
previously  Identified.  Atmospheric  pressure,  P,  Is  in  units  of  bars. 

2.5.4  RECOMMENDATIONS 

The  heat  loss  due  to  evaporation  for  the  three  models  discussed  differs  essentially  In  the 
calculation  of  the  molecular  diffusion  and  mass  transfer  at  the  interface  of  the  liquid  pool  and 
the  atmosphere.  A  recommendation  of  the  most  appropriate  evaporative  cooling  calculation 
should  then  focus  on  the  best  mass  transfer  and  dlffuslvlty  calculation. 

The  modified  Ille  and  Springer  mass  transfer  coefficient  calculation  Is  the  only  one  that 
attempts  to  estimate  the  wind  profile,  based  on  the  time  of  day,  cloud  cover,  and  terrain 
roughness.  These  parameters  strongly  influence  the  mass  transfer  coefficient,  and  thus  the 
evaporation  rate  and  evaporative  heat  loss.  Intuitively,  one  would  expect  parameters  such  as 
terrain  roughness  and  amount  of  sunlight  to  affect  the  wind  profile,  and  thus  the  mass 
transfer  rate. 

The  Schmidt  number  Is  a  fundamental  component  of  the  mass  transfer  coefficient.  The 
Schmidt  number  is  a  dimensionless  parameter  calculated  from  the  kinematic  viscosity  and 
molecular  dlffuslvlty.  Viscosity  can  be  estimated  for  pure  air,  pure  chemical  vapor,  or  a 
chemical  vapor/alr  mixture.  While  air  properties  may  be  sufficient  for  a  very  low  chemical 
vapor  pressure  and  pure  chemical  properties  sufficient  for  a  nearly  atmospheric  pressure  of 
the  vapor,  one  does  not  know  a  priori  which  situation  will  prevail.  Furthermore,  comparable 
amounts  of  air  and  chemical  vapor  may  exist  at  the  liquid  pool/atr  interface.  Therefore,  it  is 
best  to  calculate  the  kinematic  viscosity  of  the  chemical  vapor/alr  mixture.  Kinematic 
viscosity,  |A,  Is  obtained  from  the  dynamic  viscosity,  r|,  and  the  vapor  density,  pv  via  p.  ■  q/pv. 

The  dynamic  viscosity  of  air  and  some  pure  chemical  vapors  can  be  found  in  the  literature, 
usually  as  functions  of  the  temperature  of  the  vapor. 

The  properties  of  a  mixture  are  calculated  from  the  average  molar  fraction,  MF,  of  the  pure 
chemical  in  the  mixture.  In  the  modified  llle  and  Springer  model  the  average  molar  fraction 
is  approximately  1/2  times  the  ratio  of  the  chemical  saturation  vapor  pressure  to  the 
atmospheric  pressure. 

The  viscosity  of  the  chemical/air  mixture  Is  calculated  from:33 


^mlx 


MF  MWnehH-(l  -  MF)  MWAqalf 
MF  MWl/2  +  (1  -  MF)  MWA1/a 


(33) 


where  MW  ■  molecular  weight  of  the  chemical,  and 
MWA  ■  molecular  weight  of  the  air. 


33  Perty,  R.H.,  and  Chilton,  C.H.,  eds.  (1973)  Chemical  Engineers'  Handbook,  5th  edn., 
McOraw-Hill  Book  Company,  Inc.,  New  York. 
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The  vapor  density  of  the  mixture,  pv  mlx  Is  calculated  from: 


P 


v.mlx 


[MF  MW  +  MWA  (1  -  MF>]  P  . 

J _ _ _ _ _ _ _ _ _ 9101 

RT 


(34) 


where  all  quantities  are  as  defined  previously. 

Molecular  diffusivlty  is  another  fundamental  property  of  the  evaporating  liquid  pool, 
There  have  been  many  methods  proposed  for  calculating  diffusivlty  in  a  low  pressure  binary 
gas  system,  both  empirical  and  theoretical.  An  important  requirement  of  an  evaporation 
model  is  that  input  data  be  obtainable.  Reid  et  ala4  reviewed  several  calculations  and 
identified  the  method  of  Fuller  et  alao'ai>aa  as  the  one  yielding  the  smallest  average  error.  The 
method  of  Fuller  et  al  was  described  by  Eq,  (32),  The  molecular  diffusion  volumes  can  be 
obtained  by  summing  the  appropriate  atomic  diffusion  volumes  listed  on  page  586  of  Reid  et 
al.a4  The  other  parameters,  identified  previously,  are  easily  obtainable.  Therefore,  the  Fuller 
et  al  dlfitisivity  calculation  is  recommended  for  the  new  evaporation  model, 

In  summary,  the  new  evaporation  model  will  calculate  the  evaporative  cooling  via  Eq,  (22), 
The  mass  transfer  coefficient  calculation  of  file  and  Springer  will  be  adopted,  Viscosity  and 
vapor  density  of  the  chemical  vapor  mixture  and  the  molecular  diffusivlty  of  Fuller  et  al  will 
be  used  to  calculate  the  dimensionless  Schmidt  number. 


2,6  Sensible  Heat  Transfer  Due  to  Conduction  and  Turbulenoe  (Qh«) 

2.6.1  COMPARISONS 

The  rate  of  transfer  of  sensible  heat  per  unit  area  by  air  flowing  over  the  uniformly  heated 
pool  surface  is  driven  by  the  temperature  difference  between  the  air  and  the  chemical  pool: 


Oh’ki,(T.-TP)'W“a' 


(35) 


Ta  and  Tp  are  the  temperature  of  the  atmosphere  and  the  pool,  respectively,  and  kh  is  the  heat 
transfer  coefficient.  The  heat  and  mass  transfer  coefficients  are  related  by: 


24  Reid,  R.C.,  Prausnitz,  J.M.,  and  Poling,  B.E,  (1987)  The  Properties  qf  Oases  and  Liquids,  4th 
edn.,  McGraw-Hill  Book  Company,  Inc,,  New  York. 
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(36) 


k.  =  k  p  C  (Sc/Pr)0,667,  J  m~2  s'1  K  \ 

li  m  rv  p,v  v  7 


The  mass  transfer  roeffieient  (km),  vapor  density  (pv).  and  Schmidt  number  have  been  defined 
previously.  Cp,v  is  the  vapor  heat  capacity  at  constant  pressure.  Pr  is  the  Prandtl  number,  a 
dimensionless  heat  transfer  parameter  defined  as: 


Pr  « 


CiA 

TCONV 


137) 


where  Cp,v  (J  kg-1  K"1)  and  r|v  (kg  m_1  s-1)  are  the  vapor  heat  capacity  and  dynamic  viscosity, 
respectively.  TCONV  is  the  vapor  thermal  conductivity  (J  nr1  a-1  K"1)  . 

The  three  models  use  the  same  general  procedure  to  calculate  the  net  heat  transfer  due  to 
conduction  and  turbulence.  However,  their  heat  transfer  coefficients  differ,  Just  as  their  mass 
transfer  coefficients  differ.  Differences  in  mass  transfer  coefficient  and  Schmidt  number  have 
been  described  and  recommendations  made  (Section  2.5).  The  remaining  differences  are  in 
thermodynamic  properties  of  the  vapor,  namely,  heat  capacity,  density,  viscosity,  and  thermal 
conductivity. 

In  the  modified  file  arid  Springer  model,  thermodynamic  properties  are  calculated  for  the 
mixture  of  chemical  vapor  and  air  at  the  interface  of  the  atmosphere  and  the  pool.  The  ADAM 
model  uses  properties  of  the  pure  chemical  vapor  only,  while  the  Kawamura  and  MacKay 
model  uses  properties  of  the  air  only.  Furthermore,  the  ADAM  and  the  modified  llle  and 
Springer  models  calculate  properties  as  functions  of  temperature,  while  the  Kawamura  and 
MacKay  model  uses  constant  values  of  Hir  properties. 


2.6.2  RECOMMENDATION 

The  heat  transfer  due  to  conduction  and  turbulence  is  described  by  Eqs.  (36)  and  (36).  The 
three  models  discussed  differ  In  their  treatment  of  the  mass  transfer  process  and 
thermodynamic  properties.  These  topics  were  discussed  in  Section  2.5,  As  concluded  in 
Section  2.5,  one  does  not  know  a  priori  whether  air  or  pure  chemical  will  dominate  the  vapor 
above  the  liquid  pool.  For  this  reason,  the  new  evaporation  model  follows  the  modified  llle 
and  Springer  model  by  calculating  tire  thermodynamic  properties  of  the  chemical  vapor/air 
mixture. 

A  method  was  given  for  calculating  the  viscosity  of  a  mixture  in  Section  2.5.  The  thermal 
conductivity,  heat  capacity,  and  density  of  the  vapor  should  likewise  be  calculated  for  the 
chemical  vapor/alr  mixture.  A  method  for  calculating  the  vapor  density  of  a  mixture  was 
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given  in  Eq.  (34).  The  thermal  conductivity  of  the  chemical  vapor /air  mixture  la  calculated 
from  Eq.  (38),  given  in  Reference  23: 


TCONVM1X  -  MF  TCONV  MW1/3  +  (1  -  MF)  MWA1/3  TCONVA 

MF  MW1/3  +  (1  -  MF)  MWA1/3 


(38) 


where  TCONV 

TCONVA 

TCONVM1X 


thermal  conductivity  of  the  pure  chemical  vapor, 
calculated  as  a  function  of  temperature, 
thermal  conductivity  of  air  (set  equal  to  a  constant),  and 
thermal  conductivity  of  the  mixture. 


MF,  MW,  and  MWA  were  defined  previously. 

The  heat  capacity  of  a  mixture  is  calculated  by  simple  averaging: 


MF  Cp  MW  +  (1  -  MF)  Cp ^  MWA 

Cp.raix "  MFN.W  +  (1  -  MF)  MWA  “  (39) 


Cp,mix  *  beat  capacity  of  the  mixture, 

Cp,  -  heat  capacity  of  pure  chemical  vapor,  calculated  as  a 
function  of  temperature,  and 

Cp,a  «  heat  capacity  of  the  air,  given  as  a  constant  value. 


In  summary,  the  new  evaporation  model  will  calculate  the  heat  transfer  due  to  conduction 
and  turbulence  via  Eqs.  (38)  and  (36).  The  mass  transfer  coefficient  calculation  of  modified 
Ille  and  Springer  will  be  adopted.  Viscosity,  vapor  density,  thermal  conductivity,  and  heat 
capacity  will  be  calculated  for  the  chemical  vapor/air  mixture  at  the  pool/atmosphere 
interface,  The  molecular  dlffusivity  will  be  calculated  by  the  method  given  in  Eq.  (32). 


2.7  Evaporation  Rate 


The  mass  evaporation  rate  per  unit  area,  E,  is  calculated  assuming  the  vapor  at  the  pool 
surface/air  interface  behaves  as  an  ideal  gas: 


E 


k  MP  .  , 

kg  m  i 


(40) 
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km,  the  mass  transfer  coefficient  (m/s),  is  calculated  as  In  Eq.  (24)  and  is  a  function  of  pool 
temperature,  Tp.  The  saturation  vapor  pressure  of  the  chemical,  P,at  (Pa),  is  also  a  function  of 
the  pool  temperature  via  the  Antoine  equation:18 


log 


_BL 


p..,-a-c+t 


(41) 


where  A,  B,  and  C  are  constants  unique  to  each  chemical. 

To  get  the  overall  evaporation  rate  In  kg/s,  Eq.  (40)  Is  multiplied  by  the  surface  area,  A,  of 
the  pool.  In  the  modified  file  and  Springer  and  the  Kawamura  and  MacKay  models,  the  area  is 
input  by  the  user.  In  the  ADAM  model  it  is  calculated  from  the  volume  of  chemical  spilled, 
assuming  a  pool  depth  of  1  cm. 


2,6  Calculation  of  Liquid  Chemical  Pool  Temperature 

Many  of  the  energy  terms  comprising  Eq.  (1),  the  steady  state  energy  balance,  can  be 
expressed  as  functions  of  the  pool  temperature.  The  evaporation  models  solve  for  the  pool 
temperature  using  the  Newton-Raphson  iterative  procedure. 


3.  SUMMARY  OF  THE  NEW  EVAPORATION  MODEL 


In  thf  previous  section,  three  steady  state  energy  balance  evaporation  models  for  a  pool  of 
spilled  liquid  chemical  were  compared  and  evaluated.  A  new  evaporation  model  was  proposed 
based  on  these  evaluations.  In  this  section,  a  summary  of  the  new  evaporation  model  is 
presented.  The  new  model  is  a  composite  of  the  three  models  reviewed  in  Section  2. 

Table  3  shows  the  input  data  required  for  each  model. 
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Table  3.  Input  Data  Required  for  each  Model. 


Physical  Data 

l&S 

ADAM 

K&M 

New 

atitude,  longitude 
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emperature 


umber 
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These  data  are  supplied  by  separate  chemical  data  files  and  by  user  Input  during  program 
execution.  The  modllled  Ule  and  Springer  model  and  the  new  model  require  the  most 
chemical  data.  Unavailability  of  chemical  data  may  limit  the  use  of  these  models. 

All  of  the  models  reviewed  and  the  one  proposed  are  steady  state  energy  balance 
evaporation  models.  A  steady  state  energy  balance  evaporation  model  is  one  m  which  the  sum 
of  all  sources  of  energy  (heat)  transported  into  the  pool  exactfy  balances  the  sum  of  energy 
sources  transported  out  of  the  pool.  Many  of  these  sources  of  energy  can  be  expressed  In  terms 
of  the  liquid  chemical  pool  temperature,  which  is  calculated  iteratively.  That  calculated 
temperature  is  used  to  calculate  the  evaporation  rate. 

The  steady  state  energy  balance  Is  expressed  as: 


Sol 


+  Q*»  +  0pd  +  O 


he  +  ®«v  *  ®  grd  1 


0. 


(42) 


Terms  of  this  equation  were  defined  in  Section  1  and  described  in  detail  in  Section  2. 
Calculations  for  the  new  model  are  summarised  below, 


3.1  Heat  Due  to  Net  Solar  Radiation 

The  net  heat  due  to  solar  radiation  reaching  the  spilled  liquid  depends  on  the  angle  of  the 
sun  with  respect  to  the  pool  and  the  amount  and  type  of  cloud  cover.  The  solar  altitude  angle, 
SA,  Is  calculated  using  an  equation  from  Reference  5: 


sin  SA  ■  sin  LA  sin  D  +  cos  LA  cos  D  cos  SHA. 


(43) 


LA  is  the  geographical  latitude  of  the  pool,  D  is  the  solar  declination,  and  SHA  is  the  solar 
hour  angle, 

The  heat  from  the  net  solar  radiation  is  calculated  as  In  the  ADAM  model: 


Q„„i  ■  RS  (1  -  alb),  W/m8.  (44) 


The  liquid  pool  albedo,  alb,  is  included  to  correct  for  solar  radiation  reflected  away  from  the 
pool.  The  net  radiation  per  unit  area,  RS,  is  calculated  from  the  following:8'9'10 
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RS  -  MOrtn(SA)-M(i  _  (1  }  3.4) 

RATI  v 


(45) 


where  SA  is  the  solar  altitude  angle  and  RATI  is  the  diffuse  sky  radiation.  CT  is  the  cloud 
transmissivity  (related  to  thickness)  and  CF  is  the  fraction  of  sky  covered  by  clouds.  The 
diffuse  sky  radiation,  RATI,  is  a  function  of  solar  altitude  angle  and  is  calculated  from: 


RATI  «  0,694  +  0.00349  SA  if  19.4  s  SA  <  42,  or 

RATI  -  0.49  +  0.014  SA  if  SA  <  19.4,  or 

RATI  -  0.84  if  SA  2  42. 


(46a) 

(46b) 

(46c) 


3.2  Long  Wave  Radiation  Emitted  by  the  Atmosphere  and  the  Pool 

The  long  wavelength  radiation  emitted  by  the  atmosphere,  Q*tm.  and  the  liquid  pool,  Qpoi, 
is  calculated  from  the  Stefan-Boltzmann  law: 


Q,  -Q  .  -  e  oT4 - e  cT4, W/ma, 

"•tin  "pol  a  a  p  p 


(47) 


The  emissivity  of  the  atmosphere  and  the  pool,  e*  and  ep,  are  set  equal  to  0.75  and  0.95. 
respectively.  The  Stefan-Boltzmann  constant  a  is  equal  to  5.67  x  10"®  W  m-a  K*4.  T,  and  Tp 
are  the  air  and  pool  temperatures,  The  former  is  input  by  the  user  while  the  latter  is 
calculated  iteratively.  Qpoi  is  negative  because  heat  is  leaving  the  chemical  pool. 


3.3  Heat  Conducted  from  the  Ground 

Heat  is  transferred  from  the  ground  to  the  liquid  at  the  ground  surface.  There  is  also  a 
thermal  resistance  within  the  pool  of  liquid  in  transferring  heat  from  the  bottom  of  the  liquid 
layer  to  the  top,  where  evaporation  occurs.  In  the  new  model,  heat  conduction  through  the 
ground  and  heat  transfer  through  the  liquid  layer  are  combined  in  an  overall  ground  heat 
transfer  coefficient; 
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(48) 


U _ ,  ■ 


**  [(l/hg)  +  (I/hi)] 


Wm"1  K“\ 


where  hg  Is  the  coefficient  of  heat  conduction  through  the  ground,  and  hi  la  the  liquid  heat 
transfer  coefficient. 

The  heat  tranafer  coefficient  through  the  ground  la  calculated  using  an  equation  from 
Reference  13: 


h 


« 


J  m"a  a"1 


(49) 


where  kg  *  thermal  conductivity  of  the  ground  (J  ra"1  a*-1  R"1). 

otg  ■  thermal  diffuaivity  of  the  ground  (m3  a"1),  and  t  ■  time  after  spill  (a), 

The  user  la  asked  to  choose  the  type  of  soil  or  surface  best  describing  that  on  which  the 
chemical  has  spilled.  The  choices  are  diy  soil,  wet  soil,  diy  sand,  wet  sand,  or  concrete.  The 
model  Includes  values  for  thermal  conductivity  and  thermal  diffuaivity  as  given  in  Table  2  for 
each  of  the  ground  types,  The  time  parameter  is  Input  by  the  user. 

The  liquid  heat  transfer  coefficient  will  be  estimated  by  the  method  of  Kawamura  and 
MacKay:4 


h.  ■i.lWm*,K*l 
1  4>d 


(80) 


The  liquid  thermal  conductivity,  ki,  Is  expressed  as  a  function  of  pool  temperature.  The 
initial  depth  of  the  pool,  d,  is  input  by  the  user.  The  model  us  js  the  average  pool  depth,  which 
is  estimated  to  be  one  half  times  the  initial  pool  depth,  Kawamura  and  MacKay  have  derived 
the  following  empirical  equation4  for  the  liquid  resistance  factor  0; 


4-  1/  {1  +  exp  [-0.06  (TB  -  70)]} , 


(51) 


where  TB  is  the  chemical  boiling  point  in  degrees  C, 
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3.4  Heat  Lost  Due  to  Evaporation 

The  evaporation  rate  of  the  liquid  chemical  from  the  pool  can  be  described  by  the 
equation: 


k  M  P  .  _a  i 
E  ■  — ,  kg  m  3  •  \ 
RT 
P 


(02) 


where  k  la  the  maaa  transfer  coefficient,  M  ia  the  vapor  molecular  weight  (kg/mole),  P*»t  la  ita 
saturation  vapor  pressure  (Pa),  R  la  the  Ideal  gas  constant  (J  mol”1 1C1),  and  Tp  is  the  absolute 
temperature  of  the  pool.  The  rate  of  heat  loss  due  to  evaporation  per  unit  area  is  the  produot  of 
the  evaporation  rate  and  the  heat  of  vaporisation,  Hv: 


(03) 


The  heat  of  vaporization,  H,,  la  calculated  as  a  function  of  temperature  of  the  pool. 

The  maaa  transfer  coefficient  km  is  calculated  from  the  method  of  Sutton,1*  as  in  the 
modified  file  and  Springer  model: 


a-n  -_n 

k  -  0.0292  Sc“0,fl®7  Ua+n  X5*",  m/hr. 

ul 


(84) 


U  (m/h)  la  the  wind  speed  at  the  spill  site  and  X  (m)  la  the  pool  diamter  (width  In  downwind 
direction).  The  mass  transfer  coefficient  units  are  converted  to  m/s  by  dividing  by  3300, 

The  Schmidt  number,  Sc,  is  a  dimensionless  number  calculated  by  dividing  the  kinematic 
viscosity,  p,  of  the  air/chemical  vapor  mixture  by  the  dlfluslvity,  DV>b,  of  the  chemical  vapor 

In  the  air: 


Sc  - 


(88) 


q  is  the  dynamic  viscosity  of  the  vapor  mixture  and  pv  is  the  density  of  the  vapor  mixture. 
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The  mass  diflUslvlty  of  the  chemical  vapor  in  air  la  calculated  by  the  method  of  Fuller  at 

al.ao.ai.aa 


DV 


0.00143  T 


1.78 


•b 


• ,  cma/a. 


(56) 


The  Xy'a  are  the  molecular  diffusion  volumea.  Theae  can  be  obtained  by  summing  the 
appropriate  atomic  diffusion  volumea  Hated  on  page  586  of  Reid  et  al.as  The  other  paramatera, 
identified  previoualy,  are  eaaily  obtainable.  Atmoapherlc  preaaure.  P,  la  in  unite  of  bare. 

The  conatant  n  in  Eq,  (54)  deacrlbea  the  wind  velocity  profile.  It  la  a  function  of  the  terrain 
roughness  and  the  atmoapherlc  stability  (aolar  altitude  angle,  cloud  cover,  and  wind  apeed)  and 
la  a  very  aenaitlve  parameter.  Smedman-Hogstrom  and  Hogstrom19  have  eatlmated  n. 


3.5  tenalble  Heat  Tranfer  Due  to  Ooaduetlon  and  Turbuleaoe 

The  rate  of  transfer  of  aenaible  heat  per  unit  area  by  air  flowing  over  the  uniformly  heated 
pool  aurface  la  driven  by  the  temperature  difference  between  the  air  and  the  chemical  pool; 


k|,(T,-T),W/mJ. 


(57) 


T„  and  Tp  are  the  temperature  of  the  atmoephere  and  the  pool,  reepectlvely.  and  kh  1*  the  heat 
tranefer  coefficient.  The  heat  and  maee  transfer  coefficient*  are  related  by: 


km  Py  cp,  (8c/Pf)0'“7.  ^  •"*  K"‘. 


(58) 


The  maaa  transfer  coefficient  (km),  vapor  density  (pv),  and  Schmidt  number  have  been  defined 
previously.  Cp,v  is  the  vapor  heat  capacity  at  constant  pressure  of  the  all/ chemical  mixture.  It 
la  calculated  as  a  function  of  pool  temperature.  Pr  is  the  Prandtl  number,  a  dimensionless 
heat  transfer  parameter  defined  as  follows: 


Pr  - 


CP,V 

TCONV ' 


(50) 
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where  Cp,v  (J  kg-1  K_1)  and  qv  (kg  nr1  a-1)  are  the  vapor  heat  capacity  and  dynamic  viscosity, 
respectively,  TCONV  Is  the  vapor  thermal  conductivity  (J  m"1  a-1  K"1).  Each  of  these 
thermodynamic  parameters  is  a  function  of  the  pool  temperature. 

Note  that  when  the  pool  is  warmer  than  the  air  there  will  be  a  loss  of  heat  from  the  pool  by 
this  heat  transfer  mode,  and  the  calculated  sensible  heat  flux  into  the  pool  will  be  negative. 


3.0  Pool  temperature  Calculation 

The  sources  of  heat  going  into  and  out  of  the  pool  are  summarised  as  in  Eq.  (1).  Many  of 
these  terms  are  functions  of  pool  temperature,  The  pool  temperature  is  calculated  using  the 
Newton-Raphson  iterative  method.  One  starts  with  an  initial  root  and  proceeds  until  a  certain 
convergence  criterion  is  met.  In  the  new  model,  that  criterion  will  be  met  when  the  calculated 
temperature  is  within  0.1  K  of  the  previous  temperature,  The  initial  root  is  set  equal  to  the 
boiling  temperature  of  the  chemical  or  the  ground  temperature,  whichever  is  lower,  The 
ground  temperature  can  be  taken  to  be  equal  to  the  air  temperature  unless  other  information 
is  available. 

The  equations  of  the  new  evaporation  model  are  generally  valid  between  the  freezing  point 
and  the  boiling  point  of  the  chemical.  If  the  calculated  temperature  drops  below  the  freesing 
point  of  the  chemical,  the  pool  temperature  is  set  equal  to  the  freesing  temperature.9  The 
evaporation  rate  is  then  calculated  using  Eq.  (40),  This  will  give  a  conservative  estimate  of  the 
evaporation  rate, 

If  the  calculated  pool  temperature  is  grester  than  or  equal  to  the  boiling  point  temperature, 
the  pool  temperature  is  set  equal  to  the  boiling  point  temperature.  The  liquid  in  the  pool  will 
not  attain  a  temperature  higher  than  the  boiling  point  at  normal  atmospheric  pressures.  The 
calculated  saturation  vapor  pressure  will  be  about  1  atmosphere  at  the  boiling  point,  but  may 
not  be  exactly  1  atmosphere.  To  avoid  problems  such  as  calculating  a  chemical  mole  fraction 
greater  than  1,  an  alternate  equation  is  used  to  calculate  the  boiling  temperature  evaporation 
rate  (Eq.  (3)|: 


E  -  A ,  kg/s. 

H 


(60) 


where  Hv  is  the  heat  of  vaporisation. 
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3.7  Chemical  Data  Bate 

Chemical  data  for  the  new  model  calculations  are  stored  In  a  data  file  separate  from  the 
program.  The  six  chemical  compounds  Included  In  the  new  evaporation  model  data  base  are 
Hated  In  Table  4  along  with  their  chemical  formulas  and  any  abbreviations  used.  With  the 
exception  of  the  hydrazines,  most  of  the  physical  and  thermodynamic  data  and  calculations 
were  taken  from  the  ADAM  model  data  base.  The  saturation  vapor  pressure  Is  calculated  using 
the  Antoine  equation  for  the  non-hydrazlne  chemicals.  This  data  was  obtained  from  "Lange's 
Handbook  of  Chemistiy."1® 


Table  4.  Chemicals  in  New  Evaporation  Model  Data  Base. 


Name 

Formula 

Abbreviation 

Nitrogen  Tetroxtde 

N2O4 

M 

Hydrazine 

_ NgH* _ 

• 

Monomethylhydrazine 

MMH 

Unsymmetrlcal  Dlmethylhydrazine 

UDMH 

Sulfur  Dioxide 

SOa 

Phosgene 

COCLa 

- 

Data  and  thermodynamic  calculations  for  the  hydrazines  were  obtained  from  a  variety  of 
sources.  Data  for  molecular  weight,  boiling  and  freezing  temperatures,  liquid  thermal 
conductivity,  and  saturation  vapor  pressure  were  obtained  from  Schmidt,1 a8a  reference  devoted 
exclusively  to  hydrazine  and  its  derivatives.  Vapor  viscosity  was  estimated  by  the  method 
given  in  the  modlfled  Ule  and  Springer  model  at  298  K,  A  single  temperature  estimate  was 
made  because  of  the  dissimilarity  of  the  calculation  to  the  vapor  viscosity  calculation  in  the 
new  model  and  the  lack  of  data  for  the  hydrazines,  Data  for  the  vapor  heat  capacity  and  heat 
of  vaporization  calculations  were  likewise  taken  from  the  modified  file  and  Springer  model. 
Data  for  the  vapor  thermal  conductivity  calculation  was  found  only  for  hydrazine  in  Reid  et 
al,  34  A  constant  vapor  thermal  conductivity  of  0.015  J  m_1  s"1  K_1  is  included  in  the  data  base 
for  MMH  and  UDMH.  This  number  was  chosen  based  on  the  thermal  conductivity  of 
hydrazine  and  some  organic  compounds  at  298  K.  Users  of  the  new  model  should  be  aware 
that  this  is  only  an  estimate  to  be  used  until  such  data  becomes  available.  Thermodynamic 
equations  are  valid  over  a  limited  temperature  range,  If  a  temperature  range  was  not 


38  Schmidt,  E.W.  (1984)  Hydrcuctne  and  Its  Derivatives:  Preparation,  Properties,  and 
Applications,  John  Wiley  and  Sons,  New  York, 


available,  the  freezing  point  and  boiling  points  are  used  as  the  lower  and  upper  bounds  of  the 
range  of  valid  temperatures. 


4.  NEW  MODEL  SENSITIVITY  STUDIES 


In  the  previous  sections,  a  new  evaporation  model  was  developed  based  on  the  evaluation 
of  three  existing  evaporation  models  -  the  modified  Ule  and  Springer  model,  the  ADAM  model, 
and  the  Kawamura  and  MaoKay  model.  In  this  section,  sensitivity  studies  of  the  new  model 
are  presented  In  support  of  the  recommendations  of  the  previous  sections. 

Parameters  In  the  new  model  are  changed  one  at  a  time  from  a  set  of  default  parameters  to 
observe  the  effect  on  calculated  evaporation  rate.  For  instance,  the  solar  radiation 
calculations  from  each  of  the  other  models  will  be  substituted  into  the  new  model  to  see  how 
they  affect  the  results.  Default  parameters  are  listed  in  Table  5.  If  conditions  are  listed  as 
cloudy  or  "clouds",  this  means  that  there  is  100  percent  cloud  cover  by  low  clouds 
(transmissivity  ■  0.17).  Day  refers  to  12  noon,  and  night  refers  to  12  midnight. 


Table  0.  Default  Values  of  Parameters  for  Sensitivity  Studies. 


Chemical 

Nitrogen  Tetroxide 

Initial  Pool  Temperature 

290  K 

Air  Temperature. 

290  K 

Ground  Temperature 

290  K 

10  mWind  Speed 

4  m/s 

Month 

4 

Day 

20 

Time  -  Day 

12:00 

Time  -  Night 

Volume  Spilled 

22  m3 

Average  Depth  of  Pool 

1  cm 

Cloud  Cover  Fraction 

0  (no  clouds)  or  1  (clouds) 

Cloud  Type  (for  clouds  ■  1) 

Low  Clouds 

Type  of  Terrain 

Sparse  Woods 

Wet  Soil 
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4.1  Solar  Altitude  Angle 


In  Table  1 ,  it  was  shown  that  the  differences  between  calculated  solar  angles  for  the  ADAM 
and  the  Kawamura  and  MacKay  models  are  trivial  (less  than  0.5  degree).  A  worst  case  scenario 
is  examined  in  Table  6,  which  shows  that  differences  in  solar  altitude  angle  as  high  aa  5 
percent  make  a  difference  in  the  calculated  evaporation  rate  of  around  1  percent,  Therefore, 
the  selection  of  the  ADAM  model  solar  angle  calculation,  which  was  made  based  on  the 
convenience  of  the  calculation,  does  not  bias  the  results. 


Table  6.  Effect  of  Solar  Altitude  Angle  on  Evaporation  Rate. 


Solar  Altitude 

Angle 

Temperature 
of  Pool,  K 

Evaporation  Rate 

60 

264,4 

9.9 

63 

264.5 

10.0 

66 

264,6 

10.1 

69 

264,7 

10.2 

72 

264.8 

10.3 

4.3  Net  Heat  From  8olar  Radiation 

Comparisons  of  the  three  calculation  methods  discussed  in  Section  2.1  are  shown  in 
Table  7.  As  expected,  the  different  solar  radiation  calculations  yield  the  same  evaporation 
rates  for  nighttime.  The  ADAM  model  calculation  was  incorporated  into  the  new  evaporation 
model  because  it  was  the  only  solar  radiation  calculation  to  account  for  both  amount  and 
thickness  (transmissivity)  of  the  cloud  cover.  Table  7  shows  that  for  daytime  conditions,  both 
amount  and  type  of  cloud  cover  have  some  effect  on  the  calculated  equilibrium  pool 
temperature  and  evaporation  rate.  While  that  effect  is  quite  small,  especially  for  cloud 
thickness,  it  is  still  reasonable  to  use  the  ADAM  model  solar  heating  calculation  in  the  new 
model,  The  evaporation  rates  of  the  new  model  were  most  similar  to  those  calculated  using 
the  Kawamura  and  MacKay  solar  radiation.  While  it  looks  as  though  the  modiffed  Ille  and 
Springer  model  solar  radiation  does  not  depend  on  cloud  cover,  the  evaporation  rates  are 
identical  as  a  result  of  stopping  the  calculation  when  it  reaches  the  freezing  point. 
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Table  7.  Comparison  of  Methods  for  Calculating  Net  Heat  from 
Solar  Radiation. 


Source  of  Solar  Radiation  Calculation 

Evaporation  Rates,  kg/s 

Conditions 

ADAM 

(and  New  Model) 

Kawamura  and 

MacKay 

Modified  Ille 

and  Springer 

Day, 

No  clouds 

6,9 

10.0 

8.4 

Day, 

50%  low  clouds 

9.6 

9.2 

8.4 

Day, 

100%  low  clouds 

8.4 

8.4 

8.4 

Day, 

100%  high  clouds 

8.6 

8.4 

8.4 

Night, 

No  clouds 

2.3 

2.3 

2.3 

Night, 

100%  low  clouds 

8.3 

5,3 

5.3 
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4.3  Long  Wavelength  Radiation  Emitted  by  the  Atmosphere 


The  new  model  was  run  using  a  variety  of  input  conditions  for  atmosphere  emissivities  of 
0.75  and  0,81  (Table  8),  These  represent  a  reasonable  range  of  emissivities,  as  shown  in 
Section  2,2.  In  all  cases,  the  difference  in  emissivlty  produced  an  insignificant  (<  3  percent) 
change  in  the  calculated  evaporation  rate,  Therefore,  the  choice  of  0.75  for  the  new  model 
atmosphere  emissivlty  is  reasonable. 


Table  8.  Calculations  for  Ille  and  Springer  Atmosphere 
Emissivlty  (0.75)  vs  Kawamura  and  MacKay  Atmosphere 
Emisslvily  (0.81). 


Conditions 

Atmosphere  Emissivlty 
-0,75 

Atmosphere  Emissivlty 
-0.81  . 

Temperature 
of  Pool,  K 

Evaporation 
Rate,  kg/s 

Temperature 
of  Pool,  K 

Day, 

No  Clouds 

264,3 

9.9 

264.5 

10.0 

Night, 

No  Clouds 

276.5 

2.3 

276,8 

2.3 

Day, 

Clouds 

262.0 

8.4 

262.0 

8.4 

Night, 

Clouds 

263.5 

5.3 

263.7 

S3 
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4.4  Long  Wavelength  Radiation  Emitted  by  the  Pool 


The  emissivity  of  water  ranges  from  0.9B  to  0.963  (Reference  12).  The  llle  and  Springer, 
and  ADAM  models  use  pool  emisslvitles  of  0.95  while  the  Kawamura  and  MacKay  model  uses  a 
pool  emissivity  of  0.97,  Calculated  evaporation  rate  is  not  expected  to  vary  appreciably  over 
this  narrow  range  of  pool  emissivity  values.  As  Table  9  shows,  the  difference  in  pool 
emissivity  yields  a  negligible  change  in  the  calculated  evaporation  rate.  Therefore,  the  choice 
of  0.90  for  the  new  model  atmosphere  emissivity  is  reasonable. 


Table  9.  Calculations  for  llle  and  Springer  Pool  Emissivity 
(0,95)  vs  Kawamura  and  MacKay  Pool  Emissivity  (0.97). 


Conditions 

Pool  Emissivity 
-0.95 

Pool  Emissivity 
-0.97 

Temperature 
of  Pool,  K 

Evaporation 
Rate,  kg/s 

Temperature 
of  Pool,  K 

Evaporation 
Rate,  kg/s 

Day, 

No  Clouds 

264.3 

9.9 

264.5 

9.9 

Night, 

No  Clouds 

276.5 

2.3 

276.3 

2.2 

Day, 

Clouds 

262.0 

8.4 

262.0 

8.4 

Night, 

Clouds 

263.5 

5.3 

263.4 

5.3 
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4.B  He»t  Conducted  from  the  Qround 

A  feature  of  the  new  evaporation  model  is  an  option  to  choose  the  ground  or  spill  surface 
material.  Previous  models  assumed  a  ground /surface  material  and  associated  thermal 
properties.  Table  10  shows  the  effect  on  the  pool  temperature  and  evaporation  rate  of  varying 
the  Input  ground  material.  During  the  daytime,  ground  material  makes  little  difference 
because  other  energy  exchanges  dominate.  However,  during  the  night,  heat  conduction  from 
the  ground  becomes  more  Important,  so  the  choice  of  ground  material  has  a  significant  effect 
on  calculated  evaporation  rate.  For  the  default  nighttime  conditions  considered  here, 
evaporation  rate  ranged  from  0.9  to  2.3  kg/s  for  a  2200  square  meter  spill. 


Table  10.  Effect  of  Qround  lype  on  Evaporation  Rate. 


Daytime  - 

Nighttime  - 

Qround 

No  Clouds 

No  Clouds 

type 

Temperature 

Evaporation 

Temperature 

Evaporation 

of  Pool,  K 

Rate,  kg/s 

of  Pool,  K 

Rate,  kg/s 

Dry  Soil 

202.0 

6.4 

206.3 

1.2 

Wet  Soil 

204.3 

9.9 

270.5 

2.3 

Dry  Sand 

262.0 

6.4 

203,9 

0.9 

Wet  Sand 

202.0 

6.4 

269.0 

1.4 

Concrete 

262.0 

8.4 

270.4 

1.4 
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In  Section  2,4  a  discrepancy  between  the  Kawamura  and  MacKay  and  the  ADAM  heat 
conduction  coefficient  was  presented.  It  was  concluded  that  heat  conduction  from  the  ground 
is  best  estimated  from  the  equation  for  heat  conduction  In  a  semi-infinite  solid.13  as  In  the 
ADAM  model.  In  Table  1 1 ,  a  comparison  of  the  new  model  and  the  Kawamura  and  MacKay 
model  shows  that  at  night  the  heat  conduction  from  the  ground  Is  significant  enough  that 
halving  the  value  for  the  heat  conduction  coefficient  can  decrease  the  calculated  evaporation 
rate  by  as  much  as  1/3. 


Table  1 1 .  ADAM  Ground  Heat  Transfer  Coefficient  (ADAM  h^) 
vs  Kawamura  and  MacKay  Ground  Heat  Transfer  Coefficient  iKM  hgrd), 


Conditions 


Day, 

No  Clouds 


Night, 
No  Clouds 


Day, 

Clouds 


Night, 

Clouds 


ADAM  hgrd 


Temperature 
of  Pool,  K 


264.3 


276.5 


262.0 


KM  hgrd 


263.0 


Temperature 
of  Pool,  K 


267.5 


280.2 


264.6 


268.4 


In  Section  2.4  it  was  shown  that  thermal  resistance  of  the  liquid  as  calculated  by  the 
Kawamura  and  MacKay  method  is  non-negliglble  compared  with  the  heat  conduction  through 
a  semi-infinite  solid  (that  Is,  the  ground).  At  an  average  pool  depth  of  0.0  cm,  the  liquid 
thermal  resistance  lowers  the  overall  heat  conduction  coefficient  by  7  percent.  As  shown  in 
Table  12,  this  changes  the  evaporation  rate  no  more  than  3  percent.  However,  for  a  2.5  cm 
average  pool  depth  (0  cm  initial  depth),  the  liquid  thermal  resistance  decreases  the  calculated 
evaporation  rate  by  as  much  as  14  percent,  The  only  additional  Information  required  to 
calculate  the  liquid  thermal  resistance  is  the  chemical  boiling  temperature  end  the  liquid 
thermal  conductivity.  There  are  methods  available  for  estimating  liquid  thermal 
conductivity,  usually  from  reduced  temperature.34 
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Table  12.  Thermal  Conductivity  of  the  Ground:  Ground  and 
liquid  Thermal  Resistance  vs  Ground  Thermal  Resistance  Only. 


Conditions 
(Initial  Pool 
Depth) 

hi  +  hg 

hg  Only 

Temperature 
of  Pool,  K 

Evaporation 
Rate,  ka/s 

Temperature 
of  Pool,  K 

Evaporation 
Rate,  kg/s 

WmSSWM 

264.3 

9,9 

264.3 

10.1 

Night. 

No  Clouds 
(1  cm) 

278.6 

2.3 

276.0 

2.3 

262.0 

8.4 

262.0 

8.4 

II 

■mm 

263.6 

6.3 

264.9 

5.5 

Day, 

No  Clouds 
(0  cm) 

262.3 

1.9 

233.5 

2.1 

273.9 

0,46 

278.6 

0.52 

■ 

262.0 

1.9 

262.0 

1.9 

Night. 
Clouds 
(0  cm) 

262.0 

1.1 

262.9 

1.2 
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4.6  Heat  Lose  Due  to  Evaporation 


4.6. 1  PROPERTIES  OF  VAPOR  ABOVE  POOL 

The  new  evaporation  model  uses  properties  of  the  chemical  vapor /air  mixture  rather  than 
properties  of  the  pure  chemical  vapor  (as  In  the  ADAM  model)  or  pure  air  (as  in  the  Kawamura 
and  MacKay  model).  The  consequences  of  this  decision  are  shown  in  Table  13.  Evaporation 
rates  calculated  from  using  pure  chemical  properties  can  be  more  than  twice  the  evaporation 
rates  of  the  chemical/air  mixture.  The  evaporation  rates  calculated  from  assuming  the  vapor 
above  the  pool  is  pure  air  are  somewhat  lower  than  the  evaporation  rates  for  the  mixture.  As 
one  does  not  know  a  priori  whether  air  or  chemical  will  dominate  at  equilibrium,  It  is 
Important  to  use  properties  of  the  chemical/air  mixture  In  the  calculations. 


Table  13.  Evaporation  Rates  Calculated  for  Chemical/Air  Mixture 
vs  Air  Only  vs  Pure  Chemical  Only. 


Conditions 

Evaporation  Rates,  kg/s 

Chemlcal/Alr 

Mixture 

Pure  Air 

Pure  Chemical 

Day, 

No  Clouds 

9.9 

8.7 

17.4 

Night, 

No  Clouds 

3.3 

1.6 

3.0 

Day, 

Clouds 

8.4 

6.4 

17.4 

Night, 

Clouds 

5.3 

4.5 

9.9 
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In  the  new  model,  the  chemical  vapor  mole  fraction  la  computed  by  dividing  the  aaturatlon 
vapor  pressure  at  each  iteration  by  the  atmospheric  pressure.  The  IUe  and  Springer  model 
computes  an  average  mole  fraction  by  dividing  this  value  by  2.  Table  14  shows  the  sensitivity 
of  the  mole  fraction  calculation.  Dividing  the  instantaneous  value  by  2  has  the  effect  of 
reducing  the  calculated  evaporation  rate  by  6  percent  to  13  percent, 


Table  14,  Sensitivity  of  Pure  Chemical  Mole  Fraction  (MF) 
Calculation. 


Conditions 

MF-iaal- 

P 

atm 

ap.<» 

Temperature 
of  Pool,  K 

Evaporation 
Rate,  kg/s 

Temperature 
of  Pool,  K 

Day, 

No  Clouds 

294.3 

9,9 

295.7 

9,4 

Night, 

No  Clouds 

279,5 

2,3 

277.9 

2,0 

Day, 

Clouds 

292.0 

6,4 

292,0 

7.3 

Night, 

Clouds 

293.0 

0.3 

294.9 

0.0 

4.0.2  MASS  TRANSFER  COEFFICIENT 

The  options  for  calculating  the  mass  transfer  coefficient  were  compared  in  Section  2.5. 

The  modified  Ule  and  Springer  mass  transfer  coefficient  was  incorporated  into  the  new  model. 
It  accounts  for  the  effect  of  terrain  type  and  atmospheric  stability  on  the  wind  profile  near  the 
ground  through  the  use  of  a  stability  parameter  n.  Table  15  shows  the  aenaitivlty  of  the 
stability  parameter.  The  calculated  evaporation  rate  ranges  from  3  to  86  kg/a  for  the  default 
conditions  over  the  range  of  possible  stability  parameter  values.  By  contrast,  the  Kawamura 
and  MacKay  model  assumes  a  constant  n  value  of  0.25,  which  corresponds  to  an  evaporation 
rate  of  17  kg/s  for  the  default  conditions. 


Table  15.  Sensitivity  of  Stability  Parameter  n. 


n 

Temperature 
of  Pool,  K 

Evaporation 

Rate,  kg/s 

Mass  Transfer 
Coefficient,  m/s 

0.077 

262* 

86 

0,0494 

0.095 

262* 

71 

0,0407 

0.182 

262* 

29 

0,0110 

0.230 

262* 

19 

0.0107 

0.319 

264.3 

9.9 

0.00800 

0.333 

265.5 

9.6 

0.00453 

0.374 

268.6 

8.5 

0.00343 

0.387 

269.7 

8.3 

0.00319 

0.400 

270.7 

8.0 

0.00289 

0.450 

274.5 

7.1 

0.00213 

0.630 

286.3 

4.7 

0.00080 

0.675 

288.8 

4.2 

0.00064 

0.701 

290.3 

4.0 

0.00066 

0.750 

292.7 

3.6 

0.00045 

0.797 

294.4 

3.2 

0.00035 

1.024 

294.4 

3.2 

0.00046 

0.25 

282 

17 

0.00978 

*  freezing  temperature 


Table  10  show*  evaporation  rate*  calculated  using  the  mass  transfer  coefficient  of  the 
modified  file  and  Springer  model  compared  with  those  calculated  using  the  ADAM  model  mass 
transfer  coefficient ,  The  stability  parameter  allows  calculation  of  evaporation  rates  that 
differentiate  between  conditions  of  day  vs  night,  cloudy  vs  no  clouds,  and  thick  va  sparse 
vegetation.  Hie  ADAM  model  cannot  distinguish  between  some  contrasting  conditions. 


Table  16.  Effect  of  Mass  Transfer  Coefficient  Calculation  on 
Calculated  Evaporation  Rate, 


Mass  Transfer  Rate  Coefficient  From: 

Conditions 

Modified  file  and  Springer  Model 
(and  New  Model) 

ADAM  Model 

Mass  Transfer 

Coefficient,  m/e 

■EMCTMli 

M  >rrT7TTPPT» 

■RMH 

o.oosoa 

9,9 

0.00498 

9.8 

Night, 

No  Clouds, 
Scares  Woods 

0,00061 

2.3 

0,00483 

8,4 

Bay, 

Clouds, 
Sparse  Woods 

0.00493 

8.4 

0.00483 

8.4 

0,00381 

0,3 

0.00483 

8.4 

Day, 

No  Clouds, 
Forest 

0.00344 

8,0 

0.00498 

9.8 

Day, 

No  Clouds, 

Crass 

0.0407 

71 

9.8 

Day. 

No  Clouds, 
Desert 

0.0495 

80 

mm 

9.8 

Table  17  shows  the  effect  of  terrain  type  on  the  evaporation  rate  through  the  stability 
parameter.  Evaporation  rate  is  very  sensitive  to  the  choice  of  terrain  type.  Evaporation  rates 
for  desert  terrain  are  4  to  14  times  higher  than  evaporation  rates  for  forest  terrain. 
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Table  17.  Effect  of  Terrain  Type  on  Calculated  Evaporation  Rate, 


Conditions 

Type  of  Terrain 

Forest 

Flat/Qrass 

Desert 

Day, 

No  Clouds 

8,5  kg/s 

0.9  kg/s 

71  kg/s 

86  kg/s 

Night. 

No  Clouds 

1.0 

2.3 

5.1 

7,4 

Day, 

Clouds 

6,0 

8.4 

71 

86 

Night, 

Clouds 

4.3 

5.3 

19 

29 

Table  18  shows  the  effect  of  wind  speed  on  the  evaporation  rate,  both  directly  and  through 
the  (stability  parameter  (In  Eq,  (86)  for  the  maaa  transfer  coefficient].  The  fact  that  daytime 
evaporation  rates  for  8  m/s  winds  are  lower  than  those  for  4  m/s  winds  Indicates  that  the 
stability  parameter  (which  depends  on  the  terrain,  sunlight,  and  wind  speed)  can  actually  be 
more  important  than  the  wind  speed  alone  in  determining  the  evaporation  rate. 


Table  18.  Effect  of  Wind  Speed  on  Calculated  Evaporation  Rate. 


Conditions 


1  m/s 


Wind  Speed 


4  m/s 


8  m/s 


12  m/s 


Day, 

No  Clouds 


7.3 


9.9 


9,6 


10.4 


Night. 
No  Clouds 


0.36 


2.3 


7.6 


10.0 


Day, 

Clouds 


4,8 


8.4 


7,6 


10.0 


Night, 

Clouds 


1,6 


5.3 


7.6 


10,0 


a.  CONCLUSIONS 


A  new  evaporation  model  has  been  presented,  incorporating  elements  of  several  existing 
evaporation  models.  The  users  of  the  existing  models  need  to  know  how  different  the  results 
of  their  models  are  from  the  new  model  results.  Tables  19-21  shows  comparisons  of  calculated 
evaporation  rates  from  the  existing  models  and  the  new  model.  Default  values  in  Table  6  were 
used. 

The  Kawamura  and  MacKay  model  results  shown  in  Table  19  make  less  distinction 
between  the  day  and  night  conditions  than  does  the  new  evaporation  model.  Evaporation  rate 
is  directly  proportional  to  the  mass  transfer  coefficient  (see  Eq.  (22)].  However,  the  mass 
transfer  coefficient  in  the  Kawamura  and  MacKay  model  does  not  include  cloud  cover  or  solar 
angle.  Therefore,  the  evaporation  rates  calculated  by  Kawamura  and  MacKay  depend  more  on 
calculated  temperature  and  saturation  vapor  pressure.  Terrain  type  is  not  used  in  the 
Kawamura  and  MacKay  model  to  calculate  the  evaporation  rate.  As  indicated  by  Tables  17  and 
19.  the  Kawamura  and  MacKay  model  results  are  most  similar  to  new  model  results  calculated 
for  sparse  woods.  The  Kawamura  and  MacKay  nighttime  evaporation  rates  are  overestimated 
compared  with  the  new  model. 


Table  19.  Comparison  of  Kawamura  and  MacKay  Model  with  New  Model. 


■H 

Cloud 

Cover  % 

Cloud  Terrain 

Type  type 

(New  Model  Only) 

Evap.  Rate, 

kg/9,  K&M 

Evap.  Rate, 
kg/s,  New 

Day 

0 

n/a 

Sparse 

Woods 

16.1 

9.9 

Night 

0 

n/a 

Sparse 

Woods 

13.2 

2.3 

Day 

100 

Low 

Sparse 

Woods 

14.3 

8.4 

Night 

100 

Low 

Sparse 

Woods 

13.5 

5.3 

Day 

50 

Low 

Sparse 

Woods 

15.7 

9.6 

Day 

100 

High 

Sparse 

Woods 

14.3 

8.6 

The  ADAM  evaporation  model  (Just  *  email  part  of  the  ADAM  model)  yields  evaporation 
rates  approximately  2  to  10  times  higher  than  the  new  evaporation  model  for  conditions 
shown  in  Table  20.  Furthermore,  ADAM  evaporation  rates  vary  little  from  day  to  night  and 
from  cloudy  to  cloudless  conditions.  Terrain  type  Is  not  used  in  calculating  ADAM 
evaporation  rates. 

Table  20.  Comparison  of  ADAM  Evaporation  Model  with  New  Model. 


Table  21.  Comparison  of  Modified  Hie  and  Springer  Model 
with  New  Model. 


Time 
of  Day 

Cloud 

Cover 

% 

Cloud 

Type 

(1% 

Cloud 

<8T> 

n 

Celling. 

m 

(l&S) 

Terrain 

Type 

Evap. 

Rate. 

kg/s. 

(I*S) 

icvap. 

Rate, 

kg/s, 

(New) 

Day 

0 

n/a 

n/a 

n/a 

Sparse 

Woods 

17 

9.9 

Night 

0 

n/a 

n/a 

n/a 

sparge 

Woods 

2.4 

2.3 

Day 

100 

Stratus 

Low 

1000 

sparse 

Woods 

8.4 

8.4 

Day 

100 

Stratus 

Low 

200 

Sparse 

woods 

8.4 

8.4 

Day 

100 

Cirrus 

High 

12000 

Sparse 

Woods 

14 

8.6 

Day 

100 

Cirrus 

High 

6000 

Sparse 

Woods 

14 

8.6 

Night 

100 

n/a 

Low 

n/a 

sparse 

Woods 

7.2 

B.3 

Night 

100 

n/a 

High 

n/a 

Sparse 

Woods 

7.2 

B.3 

Day 

90 

n/a 

Low 

1000 

Sparse 

Woods 

14 

8.4 

Day 

90 

n/a 

Low 

200 

Sparse 

Woods 

14 

8.4 

Day 

90 

n/a 

High 

12000 

Sparse 

Woods 

17 

9.0 

Day 

90 

n/a 

High 

8000 

Sparse 

Woods 

17 

9.0 

Day 

BO 

n/a 

Low 

n/a 

Sparse 

Woods 

17 

9.6 

Day 

BO 

n/a 

High 

n/a 

Sparse 

Woods 

17 

9.8 

Day 

0 

n/a 

n/a 

n/a 

Desert 

BO 

86 

Night 

0 

n/a 

n/a 

n/a 

Desert 

9.1 

7.4 

Day 

100 

Stratus 

Low 

1000 

Desert 

21 

86 

Day 

100 

Stratus 

Low 

200 

Desert 

21 

66 

Night 

100 

n/a 

Low 

_ 2Z» _ 

Desert 

20 

30 

40 


Comparisons  of  the  new  model  with  the  Hie  and  Springer  model  are  shown  in  Table  21. 

The  new  evaporation  model  is  most  similar  to  the  modified  Ule  and  Springer  model.  These 
are  the  only  evaporation  models  presented  which  use  the  stability  parameter  to  estimate  the 
effect  of  the  wind  profile.  A  drawback  of  the  modified  Ille  and  Springer  model  is  that  when 
cloud  cover  is  less  than  100  percent  and  the  ceiling  is  greater  than  4878  m,  or  when  the  cloud 
cover  is  &  50  percent,  the  model  calculates  the  same  evaporation  rate  as  for  no  cloud  cover. 

The  modified  Ille  and  Springer  model  does  not  consider  cloud  type  when  cloud  cover  is  less 
than  100  percent,  but  it  does  allow  selection  of  a  ceiling  height,  which  is  related  to  cloud  type. 
It  also  considers  terrain  type,  as  indicated  by  the  difference  between  sparse  woods  and  desert 
terrain.  The  modified  Ille  and  Springer  model  yields  higher  evaporation  rates  than  the  new 
model  for  the  sparse  woods  terrain.  It  yields  mostly  lower  evaporation  rates  than  the  new 
model  for  desert  terrain.  In  general,  the  new  model  yields  a  greater  distinction  between  cloud 
conditions  and  terrain  types, 

The  new  evaporation  model  presented  in  this  report  is  a  composite  of  existing  steady  state 
evaporation  models,  and  is  an  improvement  upon  existing  models,  The  model  is  simple 
enough  to  run  on  a  PC,  but  is  complex  enough  to  require  a  substantial  amount  of  chemical 
data. 

To  put  the  model  results  into  perspective,  it  is  worth  noting  some  of  its  shortcomings,  The 
new  model  assumes  that  the  spill  surface  is  fiat  and  impermeable.  That  is  unlikely  to  be  the 
case,  unless  the  spill  surface  is  concrete,  The  spill  is  assumed  to  be  instantaneous,  which  is 
not  always  true,  The  spill  surface  area  is  assumed  to  be  constant,  and  an  average  uniform 
depth  is  assumed.  Again,  these  are  unlikely  to  hold  true  in  a  real  spill.  In  calculating  the  heat 
due  to  solar  radiation,  a  constant  pool  albedo  is  assumed.  In  reality,  the  albedo  of  the  pool 
will  vary  with  solar  angle.  In  a  spill  area  with  dense  vegetation,  buildings,  or  other  structures, 
the  pool  may  be  shielded  from  direct  sunlight,  even  on  a  cloudless  day.  Furthermore,  season 
will  determine  the  sparseness  of  the  vegetation,  This  affects  both  terrain  roughness  and  solar 
radiation  reaching  the  pool.  The  stability  parameters  used  in  calculating  the  mass  transfer 
coefficient  are  from  a  single  reference  source,  and  apply  only  to  rural  areas.  An  independent 
verification  of  these  data  would  place  greater  confidence  in  them. 

The  new  model  does  not  consider  temperature  changes  once  the  temperature  is  greater  than 
the  boiling  point  or  less  than  the  freezing  point.  It  is  unlikely  that  the  liquid  chemical 
temperature  will  exceed  the  boiling  point  at  ambient  pressures.  However,  the  temperature 
could  continue  to  decrease  once  the  pool  has  frozen.  A  number  of  factors  are  affected  if  the 
chemical  freezes  and  falls  below  the  freezing  point.  The  heat  of  sublimation  (rather  than  the 
heat  of  vaporization)  would  be  required  to  calculate  the  evaporation  rate.  The  ground  thermal 
properties  may  change  if  the  ground  is  also  frozen.  The  thermal  resistance  of  the  pool  itself 
will  change  when  it  becomes  a  solid.  Both  of  these  affect  the  heat  conducted  from  the  ground. 
The  albedo  of  a  frozen  pool  is  much  higher  than  that  of  a  liquid  pool,  so  the  amount  of  direct 
solar  heating  would  be  less  than  for  a  liquid.  The  vapor  pressure  equation  is  generally  not 
valid  for  a  solid.  A  temperature-vapor  pressure  relationship  for  the  solid  chemical  would  be 
required, 
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